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Universal scaling behavior at the upper critical dimension 
of non-equilibrium continuous phase transitions 
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In this work we analyze the universal scaling lunctions and the critical exponents at the upper 
critical dimension of a continuous phase transition. The consideration of the universal scaling 
behavior yields a decisive check of the value of the upper critical dimension. We apply our method 
to a non-equilibrium continuous phase transition. But focusing on the equation of state of the phase 
transition it is easy to extend our analysis to all equilibrium and non-equilibrium phase transitions 
observed numerically or experimentally. 
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One of the most impressive features of continuous 
phase transitions is the concept of universahty that al- 
lows to group the great variety of different types of 
critical phenomena into a small number of universality 
classes (see Q for a recent review). All systems belong- 
ing to a given universahty class have the same critical 
exponents and the corresponding scaling functions (equa- 
tion of state, correlation functions, etc.) become identical 
near the critical point. Classical examples of such uni- 
versal behavior are for instance the coexistence curve of 
liquid-vapor systems and the equation of state in fer- 
romagnetic systems (see for instance 0, Checking 
the universality class it is often a more exacting test to 
consider scaling functions and amplitude combinations 
(which are just particular values of the scaling functions) 
rather than the values of the critical exponents. While for 
the latter ones the variations between different universal 
classes are often small the amplitude combinations and 
therefore the scaling functions may differ significantly 
(see 0). A foundation for the understanding of the con- 
cept of universality as well as a tool to estimate the val- 
ues of the critical exponents was provided by Wilson's 
renormalization group (RG) approach |^ which maps 
the critical point onto a fixed point of a certain transfor- 
mation of the system's Hamiltonian, Langevin equation, 
etc. 

Furthermore the RG explains the existence of an upper 
critical dimension above which the mean-field theory 
applies whereas it fails below Dc- At the upper critical 
dimension the RG equations yield mean-field exponents 
with logarithmic corrections [3- These logarithmic cor- 
rections make the data analysis quite difficult and thus 
most investigations are focused on the determination of 
the correction exponents (see Eqs. (|5I6() below) only, lack- 
ing the determination of the scaling functions. 

In this work we investigate the universal scahng behav- 
ior of a continuous phase transition at Dc and develop a 
method of analysis that allows us to determine the expo- 



nents as well as the scaling functions. Therefore we con- 
sider three different non-equilibrium systems exhibiting a 
continuous phase transition into an absorbing phase. Fo- 
cusing on the equation of state our method can be easily 
applied to all equilibrium as well as non-equilibrium con- 
tinuous phase transitions observed in numerical simula- 
tions or experiments (as long as the conjugated field can 
be physically realized). In all three models the dynamics 
obey particle conservation and according to the univer- 
sality hypothesis of Q all models are expected to belong 
to the universality class of absorbing phase transitions 
with a conserved field. 

The first considered model is the conserved lattice gas 
(CLG) which was introduced in In the CLG lattice 
sites may be empty or occupied by one particle. In order 
to mimic a repulsive interaction a given particle is consid- 
ered as active if at least one of its neighboring sites on the 
lattice is occupied by another particle. If all neighboring 
sites are empty the particle remains inactive. Active par- 
ticles are moved in the next update step to one of their 
empty nearest neighbor sites, selected at random. 

The second model is the so-called conserved transfer 
threshold process (CTTP) Here, lattice sites may be 
empty, occupied by one particle, or occupied by two par- 
ticles. Empty and single occupied sites are considered 
as inactive whereas double occupied lattice sites are con- 
sidered as active. In the latter case one tries to transfer 
both particles of a given active site to randomly chosen 
empty or single occupied nearest neighbor sites. 

The third model is a modified version of the Manna 
sandpile model 0| the so-called fixed-energy Manna 
model In contrast to the CTTP the Manna model 
allows unlimited particle occupation of lattice sites. All 
lattice sites which are occupied by at least two particles 
are considered as active and all particles are moved to 
the neighboring sites selected at random. 

In our simulations (see 0, 0| for details) we start 
from a random distribution of particles and all models 
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reach after a transient regime a steady state which is 
characterized by the density of active sites p^. The den- 
sity is the order parameter and the particle density p 
is the control parameter of the absorbing phase transi- 
tion, i.e., the order parameter vanishes at the critical 
density pc according to p^ oc 5p^ , with the reduced con- 
trol parameter 6p = p/ pc — I. Additionally to the order 
parameter we consider its fluctuations Ap^. Approaching 
the transition point from above {Sp > 0) the fluctuations 
diverge according to Ap^ esc 6p~'^ (see Below 
the critical density (in the absorbing state) the order pa- 
rameter as well as its fluctuations are zero in the steady 
state. 

Similar to equilibrium phase transitions it is possible 
in the case of absorbing phase transitions to apply an 
external field h which is conjugated to the order parame- 
ter, i.e., the field causes a spontaneous creation of active 
particles (see for instance [13 )• A realization of the exter- 
nal field for absorbing phase transitions with a conserved 
field was recently developed in where the external 
field triggers movements of inactive particles which may 
be activated in this way. At the critical density p^ the 
order parameter and its fluctuations scale as p^ oc 
and oc h~'^ , respectively. 

Before we focus our attention to the scaling behavior 
at the upper critical dimension we briefly reconsider 
the scaling behavior below and above D^- In both cases 
the order parameter obeys for all positive values of A the 
universal scaling ansatz 



(1) 



The universal scaling function R(x, y) is the same for all 
systems belonging to a given universality class whereas 
all non-universal system-dependent features (e.g. the 
lattice structure, the range of interaction, the update 
scheme, etc.) are contained in the so-called non-universal 
metric factors a^, a^, and [l^ . Using the transforma- 
tion A a~^^^\ the number of metric factors can be 
reduced to Cp — a^a'^^/^ and c,, = a^a'^'^/f^ . We will see 
that this simple reduction is not possible at the upper 
critical dimension D^. Thus instead of this transforma- 
tion we set in the following = 1 for D ^ Dc in order 
to formulate for all dimensions a unified universal scaling 
scheme. 

The universal scaling function R is normed by the 
conditions R{1,0) — R{0,1) = 1 and the non-universal 
metric factors can be determined from the amplitudes of 
p,{dp, /i = 0) - (flp SpY^ and p,{Sp = 0, /i) - {a^ hf" . 
These equations are obtained by choosing in the scaling 
ansatz Eq. a^SpX = 1 and a,^hX^ = 1, respectively. 
Furthermore, the choice a^hX^ — 1 leads to the well 
known scaling equation of the order parameter 

p^{5p,h) ^ {a^hf^'' R{a,Spia,hy^^''^l). (2) 

Thus plotting the rescaled order parameter {ahh)~^/" pa 
as a function of the rescaled control parameter 
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FIG. 1: The universal scaling function of order parameter 
and its fluctuations (lower right inset) for the CLG model, 
the CTTP as well as the Manna model for D = 3 with (3 = 
0.840 and a = 2.069. The values of the non-universal metric 
factors are listed in TableQ The upper left inset displays the 
non-universal scaling plots accordingly neglecting the non- 
universal metric factors. For all considered models the scaling 
plots contain at least four different curves corresponding to 
four different field values (see for [TlLfT^ details). 



ap5p{ahh)~^^'^ the corresponding data of all systems in 
a given universality class have to collapse onto the sin- 
gle curve i?(x, 1). This is shown in Fig. for the CLG 
model, the CTTP and the Manna model for = 3. In 
the case that metric factors are neglected one observes 
the non-universal scaling behavior where each model is 
characterized by its own scaling function (see inset of 
Fig.lII. 

Similar the order parameter fluctuations are expected 
to obey the scaling ansatz 



a^Ap^{Sp,h) = A'' D{ap5p X^a^hX'^). 



(3) 



Again the number of metric-factors can be reduced by a 
simple transformation to — Upu]/'' and = a,^a'J^ . 
But it is instructive to use the above ansatz [Eq. J^J] 
since exactly one new metric factor (a^) is introduced 
for the fluctuations and furthermore the universal func- 
tions R and D are characterized by the same metric fac- 
tors. Identical metric factors for R and D occur for in- 
stance naturally in equilibrium thermodynamics where 
both functions can be in principle derived from a single 
thermodynamic potential, e.g. the free energy. In the 
case of non-equilibrium phase transitions one can argue 
that both functions can be derived from a corresponding 
Langevin equation. Setting Z)(0, 1) = 1 the non-universal 
metric factor Oa can be determined by the amplitude of 
the divergence of Ap^ similar to the order parameter. In 
the inset of Fig. ^ we plot the rescaled fluctuations as a 
function of the rescaled order parameter, i.e., the univer- 
sal scaling function D(x^ 1). Similar to the equation of 
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FIG. 2: The universal scaling function of the order param- 
eter and its fluctuations (inset) above the upper critical di- 
mension Dc = 4 with (3 = 1 and cr = 2. The numerical data 
agree perfectly with the universal mean-field scaling functions 
R{x, 1) and D(x, 1) (thick dashed lines). 



FIG. 3; The universal scaling function at the upper critical 
dimension Dc = 4. The right insets show the order parameter 
at the critical density and for zero field, respectively. The or- 
der parameter is rescaled according to Eqs. H5I6|I . Approach- 
ing the transition point (/i — > and Sp — > 0) the data tend to 
the function f{x) = x (dotted lines) as required. 



state we get a good data collapse of the corresponding 
data. 

We consider now the scaling behavior above the up- 
per critical dimension D^- According to the renormal- 
ization group scenario the stable fix-point of the renor- 
malization equations is usually the trivial fix point with 
classical (mean-field) universal quantities. Thus, in con- 
trast to the situation below Dc the critical exponents as 
well as the universal scaling functions are independent 
of the particular value of the dimension for D > Dc- 
In most cases it is possible to derive these mean-field 
exponents and even the scaling functions exactly since 
correlations and fluctuations can be neglected above D^- 
The mean-field scaling behavior of the CLG model and 
the CTTP was considered in and agrees with that 
of directed percolation, i.e., the scaling functions are 
given by 13 Rix,y) = x/2 + [y + (a;/2)2]i/2 and 
D{x,y) = R{x,y)[y + (a;/2)^]^^/-^. One can easily show 
that (3 = 1, a = 2, and 7' = 0. The latter case cor- 
responds to a jump of the fluctuations at the critical 
point which was already observed in numerical simula- 
tions 0,0. 

In Fig.[21we plot the rescaled order parameter as well 
as the rescaled order parameter fluctuations for D = 5 
and D ~ 6. In all cases the numerical data are in a 
perfect agreement with the mean-field scaling functions 
R{x, 1) and D{x, 1), respectively. Thus we clearly get the 
upper bound for the critical dimension, namely Dc < 5. 
This is a non-trivial result since a recently performed 
phenomenological field theory predicts the too large value 

Dc^gM- 

We now address the question of the scaling behavior at 
the upper critical dimension Dc = 4. Here the scaling be- 



havior is governed by the mean-field exponents modified 
by logarithmic corrections. For instance the order param- 
eter obeys in leading order p^{dp,h = 0) cx (5p|ln(5/9|^ 
and p^{Sp ~ 0,h) (X ^/h\\nh\^, respectively. The log- 
arithmic correction exponents B and S are characteris- 
tic features of the whole universality class similar to the 
usual critical exponents. Thus it was rather surprising 
that recent numerical investigations of the CLG model 
(B = 0.24, E = 0.45) and of the CTTP (B = 0.15, 
S = 0.28) reveals different values of the logarithmic cor- 
rection exponents 0| . In the following we will develop a 
complete scaling scenario at the upper critical dimension 
which agrees which the RG conjecture, i.e., all considered 
models are characterized by the same critical exponents, 
the same logarithmic correction exponents as well as the 
same universal scaling functions. 

As argued in \vi\ we assume that the universal scaling 
ansatz of the order parameter obeys in leading order 

a^p4Sp,h) ~ X-f^ \lnXf R{a,6p XllnXl^a^h X" llnXl"). 

(4) 

Thus the order parameter at zero field {h — 0) and at 
the critical density ((5/? = 0) is given in leading order by 

a^p^{Sp,h = 0) ^ a>|lna>|B i?(l,0), (5) 
a^p^{Sp^O,h) - ^/^\h-l^/^f R{0,1) (6) 

with Ji — b + I and S = s/2 + I and where we use the 
mean- field values (3 ~ I and ct = 2, respectively. Similar 
to the case D ^ Dc we set again ^(0, 1) — R{1, 0) — 1. 

Although the universal scaling ansatz [Eqs. H4l6|l ] and 
the non-universal scaling ansatz (without metric factors) 
are asymptotically equal, they may lead to different re- 
sults for numerically available data. For instance the 
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non-universal metric factor in Eq. ISJ results in the cor- 
rection factor |1 -I- lnOp/ln5p|^ compared to the non- 
universal ansatz. This factor tends to one for 5p Q 
but in numerical simulations 5p is hardly smaller than 
10"'^ which explains why different values of B and E are 
observed numerically 

According to the ansatz Eq. the scaling behavior of 
the equation of state is given in leading order by 



a, p,{5p, h) 



R{x,l) 



(7) 



where the scaling argument is given in leading order by 
X = ap5p\/aji I \n\/a^h\" with S = 6 — s/2 = B — E. 
Similarly we use for the order parameter fluctuations the 
ansatz 



Ap^{5p,h) 

A^' llnAl'^ D{a,SpX\\nX\\a^h A"'" jlnAH. 



(8) 



Using the mean-field value 7' = and taking into ac- 
count that the order parameter fluctuations remain fi- 
nite at Dc we get the scaling function 

Ap^{Sp,h) ~ D(x,l) . The non-universal metric fac- 
tor Oa is determined by the condition I?(0, 1) — 1. 

Thus the scaling behavior of the order parameter and 
its fluctuations at the upper critical dimension is deter- 
mined by two independent exponents (B and S) and four 
non-universal metric factors {a^, Up, a^, a^). We deter- 
mine these values in our analysis by the following con- 
ditions which are applied simultaneously: first, both the 
rescaled equation of state and the rescaled order param- 
eter fluctuations have to collapse to the universal func- 
tions R{x, 1) and D{x, 1) for all considered models. Sec- 
ond, the order parameter behavior at zero field and at 
the critical density is asymptotically given by the sim- 
ple function f{x) = x if one plots [a^p^{6p,0)/ap6p]^/^ 



TABLE L The non-universal quantities for various dimen- 
sions. The uncertainty of the metric factors is less than 
5%. For greater uncertainties the corresponding data sets dis- 
play significant deviations from the presented universal scal- 
ing plots. 

Model D Pe eta fflp o.h 



CLG, 


3 


0.21791 ±0.00009 


1 


0.434 


0.391 


8.881 


CTTP 


3 


0.60489 ± 0.00002 


1 


0.384 


0.093 


24.51 


Manna 


3 


0.60018 ± 0.00004 


1 


0.311 


0.074 


32.24 


CLG 


4 


0.15705 ±0.00010 


4.307 


1.664 


8.021 


7.327 


CTTP 


4 


0.56705 ± 0.00003 


0.689 


0.269 


0.047 


17.18 


Manna 


4 


0.56451 ± 0.00007 


0.690 


0.245 


0.040 


18.82 


CLG 


5 


0.12298 ±0.00015 


1 


0.329 


0.665 


8.971 


CTTP 


5 


0.54864 ± 0.00005 


1 


0.461 


0.251 


18.73 


CTTP 


6 


0.53816 ± 0.00007 


1 


0.421 


0.218 


157.5 


Manna 


5 


0.54704 ± 0.00009 


1 


0.870 


0.225 


20.69 



vs. |lnap5p| and [a^p^(0, h) / a^hY-/^ vs. \ hi^/aji\, re- 
spectively. Applying this analysis we observed that con- 
vincing results are obtained for E = 0.35 and B = 0.20 
(see Tablein for the values of the non-universal scaling 
factors). The corresponding plots are presented in Fig.O 
In particular the data collapse of the equation of state 
is quite sensitive for variations of the exponents B and 
E. Thus the quality of the corresponding data collapse 
could be used in order to estimate the error-bars of the 
logarithmic correction exponents. We obtained in this 
way E = 0.35 ± 0.06 and B = 0.20 ± 0.05. 



In conclusion, the investigation of the universal scaling 
behavior presents reliable results of the logarithmic cor- 
rection exponents in contrast to the non-universal scal- 
ing analysis. Furthermore the universal scaling analy- 
sis allows to determine the value of just by checking 
whether the numerical or experimental data are in agree- 
ment with the usually known universal mean-field scaling 
functions. 
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